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The macroscopic characteristics of a sohd, such as its thermal, optical or transport prop- 
erties are determined by the available microscopic states above its lowest energy level. These 
slightly higher quantum states are described by elementary excitations (phonons, magnons, 
quasi-particles) and dictate the response of the system under external stimuli. The spectrum 
of these excitations, obtained typically from inelastic neutron[lJ and x-ray[2^ 3j scattering, 
is the spatial and temporal Fourier transform of the density-density correlation function 
(fluctuations) of the system, which dictates how a perturbation propagates in space and 
time[4j. As frequency-domain measurements do not generally contain phase information [5j, 
time-domain measurements of these fluctuations could yield a more direct method for in- 
vestigating the excitations of solids and their interactions both in equilibrium and far-from 
equilibrium. Here we show that the diffuse scattering of femtosecond x-ray pulses produced 
by a free electron laser (FEL) can directly measure these density-density correlations due 
to lattice vibrations in the time domain. We obtain spectroscopic information of the lattice 
excitations with unprecedented momentum- and frequency- resolution, without resolving 
the energy of the outgoing photon, by a Fourier transform of the femtosecond dynam- 
ics. Correlations are created via an acoustic analog of the dynamical Casimir effect [6^ 7J, 
where a femtosecond laser pulse slightly quenches the phonon frequencies, producing pairs 
of squeezed phonons at momenta +q and —q^ [9]. These pairs of phonons manifest as 
macroscopic, time-dependent coherences in the displacement correlations ^Oj that are then 
probed directly by x-ray scattering. Since the time-dependent correlations are preferentially 
created in regions of strong electron-phonon coupling, the time-resolved approach is natural 
as a spectroscopic tool of low energy collective excitations in solids, and their microscopic 
interactions, both in linear response and beyond. 

Spatial density fluctuations in nominally periodic media reduce the intensity of the Bragg 
diffraction peaks and consequently increase the weak diffuse scattering between these peaks, 
the details of which reflect the amplitudes and spatial frequencies of the fluctuations [llj. 
When these perturbations are due to lattice vibrations, the intensity of the diffuse scattering 
at momentum transfer Q is proportional to the density-density correlation (uqU^q)^ where 
Uq is the phonon amplitude at reduced wavevector q = Q — Kq and Kq is the closest 
reciprocal lattice vector to QjTTl [T2j. In typical experiments with long x-ray (or neutron) 
pulses, the measured diffuse scattering is proportional to the time-integrated mean squared 
displacements, and spectral information can only be obtained by analyzing the energy of the 
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FIG. 1. Femtosecond x-ray diffuse scattering (a) Static thermal diffuse scattering from (001) 
Ge in grazing incidence from 10 keV x-ray photons at the LCLS. Dashed boxes are the q-space 
regions shown in Fig. [sj (b) Calculated equilibrium pattern using a Born model of the forces. 
White lines indicate the boundaries of the Brillouin zones. Miller indices are also indicated, (c) 
Representative traces of the normalized change in scattering A/(t)//max induced by the optical 
laser as a function of (optical) pump- (x-ray) probe delay. 

scattered particles. Recent advances in FEL sources [131 [Hj provide sufficient flux and time- 
resolution to record snapshots of the scattered x rays on a time-scale short compared to the 
density fluctuations. As we show here in a single-crystal of the prototypical semiconductor 
germanium, temporal coherences in {uqU-q){t) induced by a sudden softening of the lattice 
yield the phonon dispersion without having to resolve the scattered photon energy. 

Figure [l] (a) shows a portion of the equilibrium x-ray diffuse scattering from a single- 
crystal of germanium at grazing incidence captured with an area detector (see supplemental 
material). The bright areas correspond to regions of reciprocal space with low frequency 
phonons that contribute strongly to the equilibrium diffuse scattering, mostly dominated by 
thermally populated acoustic phonons [121 113 • Figure [l] (b) shows the simulated thermal 
diffuse scattering from a Born model of the forces flU |T2l [15] (see supplemental information 
for additional details). The calculated pattern matches the measured diffuse scattering 
extremely well. The white solid contours in Fig. [l] (b) represent the boundaries of the 
Brillouin zones accessible in this geometry and we have also indicated the respective Miller 
indices. We note that in some speciflc situations [16], the phonon dispersion can be obtained 
from a flt of the calculated image [Fig. [l] (b)] to the experimental data [Fig. [l] (a)]. Instead, 
here we show that coherences in the displacement correlations can yield large sections of the 
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phonon dispersion directly from the measurement. 

Figure [l] (c) shows the evolution of the change in diffuse scattering intensity A/(t)//max 
induced by photoexcitation with a 50 fs infrared laser pulse centered at 800 nm. The 
two curves show the time traces for the two points labeled "u" and "v" in Fig. [l] (a), 
normalized by the maximum of the laser-off image. Photoexcitation induces an overall 
step-like increase in the scattering whose magnitude depends on momentum position, and 
oscillations at frequencies in the acoustic phonon range 1-3.5 THz. The sharpness of the 
initial step as well as the highest frequency observed ^3.5 THz were resolution-limited by 
the timing jitter in the pump-probe delay ~ 250 fs[17j. As we discuss next, we identify these 
oscillations as the resulting time-dependent correlations, {uqU-q){t)^ due to the formation of 
(squeezed) pairs of phonons from the sudden softening of the lattice forces. 

We begin by considering a sudden quench (softening) in the harmonic potential driven by 
excitation of electron-hole pairs by the laser pulse. In the specific case of the tetrahedrally 
bonded semiconductors, photoexcitation from the mainly bonding valence band states to 
the mainly anti-bonding conduction band is predicted to soften the transverse acoustic 
(TA) modes[T8]-[22j, which at high enough excitation density leads to the loss of crystalline 
order through an instability of the shear (TA) modes. In more general terms, the sudden 
quench of the potential is a consequence of the electron phonon coupling and thus is not 
restricted to the specifics of these materials. The evolution of a harmonic oscillator after a 
sudden quench of the frequency is related to vacuum squeezing, where the variance in the 
quantum displacements oscillates below the quantum limit for a quarter of the oscillator 
cycle at the expense of uncertainty in the momentum, as shown for photons [23l [24 j and 
phonons [HI |9]. This is intimately related with the dynamical Casimir effect [25] in which 
a sudden change in the boundary conditions of the electromagnetic field generates real 
photons out of the vacuum quantum fiuctuationsj^. More recently, an acoustic analog of 
this effect was demonstrated in a Bose-Einstein condensate where a sudden quench of the 
sound velocity was shown to produce correlated pairs of phonons 0. 

Assuming all phonons, with frequencies Qq and oscillator mass m, are in the ground state, 
a sudden change in the frequency Qq ^ f]'^ at t = leaves each mode in a squeezed state 
where the variance in the displacement evolves according to[27j 



= j:^ f(i + + (1 - cos(2n't)l , (1) 
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FIG. 2. Coherence in the density-density correlations, (a) Representative frames of the 
oscillatory component of A///max after background subtraction, (b) Time dependence of the 
subtracted data at a few reduced wavevector locations between u = [-0.09 -0.03 -0.08] and v = 
0.08 0.13 -0.27] (r.l.u.) in Fig. [l] (a). These curves have been displaced vertically for clarity. 

where f3q = Vtq/Vt'^ > 1 for a sudden softening. At finite temperature, Eq. ([T]) contains an 
additional thermal factor [7J. This expression describes the generation of correlated pairs 
of phonons at q and —q\J\' This mechanism is formally analogous to (spontaneous) para- 
metric down- conversion. For simplicity we have assumed that the equilibrium position of 
the oscillator is unchanged in the excited state, i.e. there is no coherent displacement and 
thus {u{t)) = 0[8j, but this is not a requirement. The momentum of the pump photons is 
small compared to the Brillouin zone dimensions, thus conservation of momentum dictates 
that the generated pairs of phonons have equal and opposite momenta. In our case, the 
softening occurs for all q and is expected to be particularly strong at the Brillouin zone 
boundary [2T1 [22], but due to the current time-resolution limit, the phonon frequencies at 
the zone boundaries were not resolved. 

The x-ray diffuse scattering intensity is directly proportional to the projection of the 
time-dependent correlations [Eq. ([T])] along Q^Uj. Accordingly, the intensity oscillates at 
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FIG. 3. Equal energy phonon momentum distribution. Magnitude of the time- Fourier trans- 
form at representative frequencies of the background-subtracted data, (top and bottom panels) 
zoomed view of the region of space labeled "1" and "2" in Fig. [T] respectively. 

2^2', with an amplitude proportional to 1 — In the limit of low-density excitation, the 
frequencies will approximate the equilibrium values, and thus the Fourier transform of the 
oscillatory component should give the phonon dispersion. In our case, |1 — ^ 0.05 and 
thus ^q/^q ^ 1.025, which is close to the frequency resolution limit of ^ 0.1 THz given 
by the finite time window in these data. Consistent with bond softening the mean square 
displacements (and thus the scattering) increases during the first quarter cycle, as seen in 
Fig.|l|(c). 

For better sensitivity to the oscillatory signal we filtered the slowly varying background 
from the raw data (supplemental information and movies). Figure [2] (a) shows representa- 
tive frames of the obtained oscillatory component. The plots in Fig. [2] show an enlarged 
view of the (022) Brillouin zone. The data for zone (113) shows similar results. The blue 
(red) regions in this figure represent an increase (decrease) in the intensity relative to the 
subtracted average. The fringes in g-space seen here originate from phonons with diflFerent 
frequencies across reciprocal space, which have phase coherence due to the sudden frequency 
softening. The traces in Fig. [2] (b) show some of these oscillations for a few wavevectors 
along the u — v segment in Fig. [l] (a). 

Figure |3] shows a zoomed view at selected frequencies of the Fourier transform (FT) along 
the time axis of the oscillatory component in Fig. [2] (top and bottom rows represent regions 
in boxes "1" and "2" in Figs, [l] (a), respectively). The intensity at each pixel and a given 
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FIG. 4. Extracted dispersion relation in selected directions, (a) Schematic representation 
of the constant-frequency cuts of the acoustic dispersion relation that yield the data in Fig. |3j 
The surfaces represent the two TA branches and the plane represents a const ant- frequency cut at 
2cl; = 2.5 THz. (b) and (c) acoustic dispersion along the sections shown with dashed hues in (d) 
where qi = [-0.09 -0.04 -0.07], q2 =[-0.04 -0.82 0.37], = [0.14 0.01 0.05] and q^ =[0.27 -0.95 
-0.07] (r.l.u). White lines in (b) and (c) represent the calculated acoustic dispersion. 

frequency is the magnitude of the FT of traces like those shown in Fig. [2] (b) . The bright loops 
appear at locations in momentum space where the intensity oscillates at the same frequency. 
These contours (Fig. |3]) represent constant-frequency cuts of the phonon dispersion relation 
as depicted schematically in Fig. [i] (a). The data in Fig. |3]show two bands, seen more clearly 
in the bottom row plots, which correspond to the two TA branches, with pinch points where 
the bands are degenerate along high-symmetry directions. Their intensity depends on the 
amplitude of the coherent mean squared displacements, as well as their projection along Q. 

Finally, Figs, [i] (b) and (c) show the extracted dispersion relation along the directions 
indicated by the dashed lines in Fig. |4] (d). Here the phonon frequency is = cl;/2 

according to Eq. ([T]). It is important to note that we have not relied on any theoretical 
model of the interatomic forces to extract the phonon frequencies, only a knowledge of the 



scattering geometry and the time-delay is required. The white hnes in Figs.|4](b) and (c) are 
the calculated equilibrium dispersion. Note that within our experimental sensitivity we only 
pick out the transverse acoustic phonon branches and not the longitudinal acoustic branch, 
(currently the resolution is not high enough to resolve either optical phonon branch) . This is 
to be expected as the excitation of carriers reduces the strength of the covalent bonds which 
give rise to the shear stability in the tetrahedrally bonded semiconductors [28j. Otherwise, 
the discrepancies are small and could be due to systematic errors in determining the sample 
orientation or the forces as much as changes in the excited state forces. The curvature of the 
branches is due to our particular geometry, which results in a non-planar section of reciprocal 
space and avoids the Bragg peak at g = where Q = 0. The flat spectral components at 
lower frequencies are likely due to fluctuations of the FEL that were not removed by our 
background subtraction. The sample was oriented far from the zone-center {q = 0) to avoid 
strong Bragg reflections on the detector, particularly given the large wavelength fluctuations 
of the FEL. Recently demonstrated self-seeded operation of the FEL|14j will provide better 
pulse stability with narrower bandwidth and will yield better momentum resolution, lower 
noise, and will allow sampling closer to g = 0. An improved timing diagnostic [29] will enable 
the observation of faster oscillations and thus higher frequencies. 

The induced temporal coherences in the density-density correlations observed here are a 
consequence of a sudden change in the interatomic potential and can thus be generalized 
to other excitations. These coherences span the entire Brillouin zone but will be favored in 
regions where the resultant (real or virtual) charge-density couples strongly to the phonons. 
For example it will be particularly strong in regions of enhanced electron-phonon coupling 
and could flnd broad use in the study of the coupled degrees of freedom in complex materials. 
We further stress that, far- from equilibrium the pump-probe approach gives unique access 
to the phonon excitations and their interactions in the short-lived transient state. 
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